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IMPACT DAMAGE OF COMPOSITE PLATES 


By 

K.M. Lai 1 and G. L. Goglia 2 
SUMMARY 

During the first half of the current year (September 1982 to March 
1983) the following work was done: 

1. Development of a model for coefficient of restitution. 

This model is based on idealized shear dominated theory of fiber- 
reinforced composite laminates which assumes that the fibers are in- 
extensible in fiber direction and that the fibers are incompressible in 
in z-direction. This model is strictly limited to thin plates of 
quasi-isotropic composite material. 

Predictions were compared with the test results for circular and 
rectangular plates of various sizes. The similarity between theory and 
experiment was found better than could reasonably be expected, in view 
of the very simple nature of the theory. 

A manuscript "Coefficient of restitution for low velocity trans- 
verse impact of thin graphite-epoxy laninates" (Appendix A) describes 
the model, and is enclosed herewith. 

2. The experience of impact tests in the laboratory revealed that the 
variation of force and central deflection during an impact event could 
be represented by half-sine waves, with a phase lag between force and 
deflection. An attempt was made to correlate the phase difference with 
the energy absorbed by the target. Since the energy absorbed by the 
target is equal to the energy lost by the impactor during impact, which 
can be easily calculated by incorporating the coefficient of restitu- 
tion, a model was developed to correlate the coefficient of restitution 
with the phase difference between force and deflection. 

Associate Research Professor, Department of Mechanical Engineering and 
Mechanics, Old Dominion University, Norfolk, Virginia 23508. 

2 Eminent Professor/ Chairman, Department of Mechanical Engineering and 
Mechanics, Old Dominion University, Norfolk, Virginia 23508. 



Comparison of the predictions with the experimental test results 
was very encouraging . 

A manuscript "Relationship between phase difference and coeffi- 
cient of restitution during low velocity foreign objects transverse 
damage of composite plates" (Appendix B) is enclosed. 

3. Work to evaluate the state of stress in various plies surrounding 
the point of impact is in process. Report is expected to come in about 
3 to 4 months . 
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APPENDIX A 
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COEFFICIENT OF RESTITUTION FOR LOW VELOCITY TRANSVERSE 
IMPACT OF THIN GRAPHITE-EPOXY LAMINATES 

K. M. Lai* 

NASA Langley Research Center 
Hampton, Virginia 


ABSTRACT 


This paper discusses a simple model to study low velocity transverse 
impact of thin plates made of fiber-reinforced composite material, in particular 
T300/5208 graphite-epoxy. This model predicts the coefficient of restitution, 
which is a measure of the energy absorbed by the target during an impact event. 
The model is constructed on the assumption that the plate is inextensible in 
the fiber direction and that the material is incompressible in the z-direction. 
Such a plate essentially deforms by shear, hence this model neglects bending 
deformations of the plate. The coefficient of restitution is predicted to 
increase with large interlaminar shear strength and low transverse shear modulus 
of the laminate. 

Predictions are compared with the test results of impacted circular and 
rectangular clamped plates. Experimentally measured values of the coefficient 
of restitution are found to agree with the predicted values within a reasonable 
error. 


INTRODUCTION 


During the last decade, the advanced graphite-epoxy composite materials 
have offered specific strength and modulus properties superior to those of 
conventional monolithic materials. These advanced composites have potential 


*Also, Research Associate Professor of Mechanical Engineering and Mechanics 
at Old Dominion University, Norfolk, Virginia. - 
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for use in aerospace structures, but they appear to be susceptible to trans- 
verse impact damage. The problem of determining the energy absorbed by the 
target is of interest, because this energy can reduce the integrity of the 
target. 

The energy absorbed by the target during an impact event is equal to the 
loss in kinetic energy of the impactor. Mathematically, the energy lost, I a , 
by an impactor during the impact process is given by 

I a = 1(1 - e 2 ) (1) 

where I is the maximum amount kinetic energy carried by the impactor and e 
is the coefficient of restitution. Thus, the coefficient of restitution is a 
measure of the energy lost by the impactor and is equal to the ratio of the 
velocity at the time of rebound to the maximum velocity of the impact. Because 
the coefficient of restitution is needed to compute the absorbed energy, the 
objective of this work is to develop a mathematical model of the coefficient of 
restitution of thin anisotropic plates. 

ANALYSIS 


Composite plates are laid up by laminating thin sheets of unidirectionally 
reinforced material. Macroscopically the fibers may be regarded as continuously 
distributed through the plate and the material is transversely isotropic about 
the axis normal to the plate. For graphite-epoxy composites, E^ » E m , i.e., 
the fiber is very much stiffer than the matrix. The resistance of a lamina to 
extension in the fiber direction may be several orders of magnitude greater 
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than its resistance to transverse extension or transverse shear. This property 


<_an be idealized by assuming that all laminae are inextensible in the fiber 
direction [l]. For mathematical convenience, it will also be assumed that 
the composite is incompressible in transverse direction. This assumption 
is frequently made in solid mechanics. 

Consider a plate having its middle surface lying m the plane z » 0. 

The fibers are in the planes z = constant and are arranged in families 
having multidirections. If <p is the angle between the x-axis and a particular 
family of fibers, the condition for inextensibility in the fiber direction is 
given by [2], 

3u _ /3u 3u„\ 2 

e x = cos 2 <p + + gj-j cos <p sin <f> + sin <p = 0 (2) 

Equation (2) will hold if the displacements in the x- and y-directions are 
zero, i.e.. 


(3) 

Uy-0 

For incompressible material, with equation (3), 3u /3z = 0, and so 

z 

u z = 0 j(x,y,t) (4) 

Here (a is the displacement in z-direction at time t; u , u , u are, 

x y z 

respectively, the displacement components in x-, y-, and z-directions. Thus 

for incompressible material, u z is uniform through the plate thickness. 

The equations (3) and (4) imply shear deformation of the plate, in the 

sense that plane sections remain plane but do not rotate. This idealized 

mode of deformation is shown m figure 1. 
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Then from equations (3) and (4) the velocity components, denoted by 


v , v , and v , are 
x* y* z* 


v x = 0, v y = 0, v 2 = 3io(x,y,t)/3t = v(x,y,t) (5) 

Equation of Motion 

Consider a plate which is struck transversely at the origin at time 
t = 0 by a mass M moving with speed V . The plate, of thickness h, is 
clamped such that to = 0 at its edges. The solution to the analogous problem 
for a beam [3] suggests that the region of the plate bounded by an outwardly 
propagating wave (curve C, shown in figure 1) moves as a rigid body and the 
region outside C is at rest. It will be assumed that the wave front remains 
concentric with the impact point at the center. It is assumed that the impact 
force is resisted by the dynamic transverse shear stress in the propagating 
wave C of radius r. The radius of the wave C, at any time t during the 
impact event, will vary in such a way that dynamic shear stresses on the cir- 
cumference of this wave balance the impact force at that moment. The equation 
of motion of the plate inside the wave front can be described by the following 
expression: 


(M + m )v 
r 


-ZiTrhSj 


(6a) 


where M is the projectile mass and m r is the mass of the region of the 
plate bounded by the wave front of radius r, which can be expressed by the 
following expression: 


m r =» TTr hp 


(6b) 
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By substitution of m r from equation (6b), the equation (6a) can be rewritten 
in the following form: 

(M + Trr^hp)v =* -2TrrhS d (6c) 

where is yield shear stress resultant on the boundary of the curve C. 

If shear deformations are assumed to be elastic-plastic, the 
resultant yield shear stress, S^, will be expressed by the 
following form: 

S d = S o + Gv n (7) 

Where S q is elastic shear stress, G is the tangential modulus of 
a strain hardening material (will be taken equal to transverse 
shear modulus), v is an equivalent strain on the boundary of the 
propagating wave front C, and n is strain-hardening exponent of 
the material. 

Considering the dynamic jump condition and discontinuity along the 
curve C, the second term in equation (7) takes the following form [2] for 
linear strain- hardening materials (n ■ 1) : 

GV 11 = pC t v (8) 

where C t is equal to /gTp, known as transverse wave velocity. In case of 
T300/5208 graphite-epoxy laminates, elastic yield stress will be taken equal 
to the interlaminar shear stress, S . Taking S q = and substituting the 
value Gv 11 from equation (8), the equation (6) becomes 

(M + 7rr^hp)v = -2Trrh(S^ + pC t v) (9) 

To solve equation (9) we rewrite, after algebraic simplifications, as 

-^■[(M + 7rr 2 hp)v] = -27rrhS i (10) 
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The operator d/dt can be written in the following form: 


d_ 

dt 


dr _ d_ 
dt ' dr 


C t dr’ 


where C t 


dr 

dt 


(ID 


Now, let 


A = (M + TTr 2 hp)v (12) 

and 

B = M + nr 2 hp 


Using the relations from equations (11) and (12), the equation (10) is expressed 
below: 


dA = _ £i dB 
t dr p dr 


(13) 


Integration of equation (13) yields 

C t A " " p B + A ’ (14) 

where A* is constant of integration. Putting the initial condition, that 
r =» 0 at v = V, in equation (14) gives the value of the constant of 
integration. A', 


A' - MCSj/p + VC t ) (15) 

Substitution of A' from equation (15) and A and B from equation (12) in 
equation (13) gives 


2 ^i 2 i ^ 

.(M + Trr z hp)v + ~ (H + TTr hp) = M j + C c V 


('16) 
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Setting the boundary condition that r = a when v = 0 in equation (16) 
gives the radius, a, of the transversely propagating wave front. Equa- 
tion (16) becomes 


a 


\JnhS. 


( 17 ) 


Equation (17) gives the radius of the propagating wave front at maximum 
deflection, which shall be called the radius of the impact affected zone. 

The mass of the transversely propagating region of the plate, m, during 
impact is given by 

2 

m = ira hp (18) 


Using the above concepts, it is possible to write the equation of motion 
during plate rebound. During rebound, it will be assumed that the deformations 
are elastic, and the resultant yield shear stress will be taken equal to the 
interlaminar shear^s tress of the laminate. Thus and v =* v r in 

equation (6) gives the equation of motion of the plate during rebound event. 

(M + 7Tr hp)v r = -ZTTrhS^ (19) 


where the suffix r represents the rebound event. Substituting the operator 
d/dt from equation (11) into equation (19) gives 



2-irhSi r 

C t M + TTr 2 hp 


Integration of equation (20) with respect to r gives 


( 20 ) 


v 

r 


log (M + irr 2 hp) + B ? 


( 21 ) 
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where B' is constant of integration. Setting the initial conditions at the 
start of rebound event, v r = 0 at r = a, gives the constant B', 


B' 


g 

log (M + Tra 2 hp) 


( 22 ) 


Putting the expression for B* from equation (22) 








°1 

v = 

r pC 


log 


M + ira hp 
2 

, M + Trr hp/ 


(23) 


Immediately after detachment, assuming that no mass of the plate vibrates after 
rebound, r = 0, the expression (23) becomes: 


v 


r 




(23a) 


From equations (17) and (18) it is evident that the mass of the propagating 
region is proportional to the impact velocity. The combined impact velocity 
of the target and projectile can be taken equal to M/M' times V where M' 
is the sum of the projectile mass and the mass of the Impact affected zone of 
the target. The substitution of V by VM/M' in equation (17) and 
v r = -eV in equation (23a) gives the relation for the coefficient of 
restitution, e, by equation (24): 


e 


- S i 
PV C, 


log 1+777 


M 

P C v V o\ 

M* 

s i ) 


(24) 


Now it is possible to calculate the total energy absorbed by the target, 

2 

I,, from equation (1), i.e., I = 1(1 - e ), where I is the maximum kinetic 
energy of the impactor. For a minimum amount of kinetic energy being absorbed 
by the target, a maximum amount of kinetic energy should be taken away by the 
impactor during its rebound; this is obtained by maximizing the coefficient of 
restitution. 
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Analytical Results 


Figure 2 shows the variation of the radius of the impact affected zone, 
a, with the impact velocity as calculated from equation (17). The value of a 
is increasing with increase in velocity V . The effect of change in shear 
modulus, interlaminar shear strength, density of the composite laminate are 
computed from equation (17) and is shown in figures 3-5. Decrease in transverse 
shear modulus G of the laminate decreases the radius of the impact affected 
zone; while the increase in interlaminar shear strength, mass density and 
thickness of the laminate lowers it. Increases in mass of the impactor and 
impact velocity cause increase in the radius of this zone as can be seen from 
equation (17). 

The variation of coefficient of restitution with impact velocity is 
shown in figure 6 as calculated from equation (24) . The value of coefficient 
of restitution is found to decrease with increase in impact velocity. Fig- 
ures 7-9 show how the shear modulus, interlaminar shear strength and density 
affect the relation between the impact velocity, V , and the coefficient of 
restitution, e. Decrease in shear modulus and density are found to increase 
the coefficient of restitution, but an increase in interlaminar shear strength 

decreases it. Thus, in order to have a low value of I , we should have a 

a* 

high value of interlaminar shear strength and a low value of shear modulus of 
the laminate. Also, an increase in interlaminar shear strength and decrease 
in shear modulus lowers the radius of the impact affected zone. 

Comparison with Test Results 

In order to compare the predictions of the coefficient of restitution, 
two types of impact tests were conducted: 
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(1) Free Drop Tests: 

Rectangular panels, clamped to provide a free impact area of 90 mm x 150 mm 
were impacted transversely by free drops of a one-inch diameter steel ball from 
different heights. The bolts were torqued uniformly to keep mounting conditions 
uniform. The free-drop impact setup is shown in figure 10. The impact-energy 
and rebound energy were obtained from the measurements of drop height and 
rebound height of the ball. Table I lists the impact test data from drop tests. 

(2) Instrumented Cantilever Impact Tests: 

A one- inch steel ball, with an accelerometer on top, was fixed on the 
free end of a fiberglass cantilever rod, figure 11. In these tests, circular 
clamped plates of 50 mm and 90 mm diameter were used. During an impact event, 
the electric signals from the accelerometer were first processed through a 
digital data-acquisition system, and then through a computer and placed onto 
cassette tapes to allow direct acceleration-time data plotting. Double inte- 
gration of the test data provided velocity- time and displacement-time data. 
Typical variation of acceleration, velocity and displacement with time is 
shown, in figure 12. Table II gives only the data on impact and rebound 
velocity because the analysis of impact force and deflection is not within 
the scope of this report. 

The experimental values of the coefficient of restitution from both types 
of test are plotted in figure 13. The line in figure 13 gives the prediction 
of the variation of coefficient of restitution with velocity of impact. All 
the test values are found to lie within a reasonable range of the predicted 
curve. Thus, this analysis provides a reasonable prediction for the coeffi- 
cient of restitution over the velocities of impact from 0 to 6 m/sec. 
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DISCUSSION 


The main concept of the approach to the dynamics of transverse loading 
of plates in this paper is the assumption of inextensibility in the fiber 
direction and incompressibility in the transverse direction. During loading 
of the impact event, the resultant shear stress around the circumference of 
the propogating zone is assumed to be the sum of shear yield stress and the 
plastic shear stress; whereas during unloading part of the event, the defor- 
mations are assumed to be elastic and the resultant shear force around the 
circumference of the propagating zone is the shear yield stress. 

In view of the idealizations made in the presented model, it is interesting 
to examine the qualitative output of the model. The simple theory provides us 
the kinetics of impact leading to a mathematical model to evaluate the coeffi- 
cient of restitution and the impact-affected zone around the point of impact. 
From equation (16) we see that area of the impact-affected zone is proportional 
to the impactor's momentum and square root of the transverse shear modulus, G, 
and inversely proportional to the interlaminar yield shear stress, S^. Thus 
small G and large lead to absorption of the kinetic energy of the 

impactor over a small area of the plate. Also low values of G and large 
values of (equation (22)) lead to higher values of coefficient of 

restitution. High values of coefficient of restitution lead to low absorption 
of kinetic energy. Thus in order to have the lowest value of the kinetic 
energy absorbed by the target (which means more nearly elastic impact) all 
the deformations must occur at the point of impact. Such a case corresponds 
to indentation impact over a plate fixed on rigid elastic plate; this case is 
outside the scope of the present model. 
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Comparison of experimental values of coefficient of restitution for 
transverse impacts of 50 mm and 90 mm diameter plates and 90 mm x 150 mm 
rectangular plates agreed well with the predictions from the present model. 
It can be concluded that in case of low velocity transverse impacts of thin 
laminates made of T300/5208, the coefficient of restitution is a material 
property which depends strongly upon transverse shear modulus, interlaminar 
shear strength and the density of the laminate. 
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SYMBOLS 


A',B' 


a 

c t 

E f’ C m 

G 

h 

I 

I 

a 

M 

M* 

m 

“r 

r 

S d 

s i 

Sw 

t 

» u z 

V 


v x’ v f v z 


w 


Constants of integration 

Maximum radius of the impact affected region of the plate 

Transverse shear wave velocity 

Young's modulii of fibers and matrix 

Transverse shear modulus 

Thickness of the plate 

Maximum kinetic energy of the impactor 

Kinetic energy absorbed by the target during impact process 
Mass of the impactor 

Mass of the impactor plus the moving region of the plate, M + m 
Maximum mass of the propagating region of the plate 
Mass of the propagating region of radius r 
Radius 

Dynamic shear stress 

Interlaminar strength of the plate 

Shear stress on the boundary of the transverse wave 

Time 

Displacement components in x-, y-, z-directions 
Maximum velocity of the impactor 
Velocity components m x-, y-, z-directions 
Vertical displacement at any x,y,t 
Mass density of the plate 

Instantaneous velocity during impact sequence 
Instantaneous velocity during rebound sequence 
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TABLE I 


TEST RESULTS OF FREE BALL DROP TESTS 


Impactor : 
Target : 

Impact area: 


1-inch diameter steel ball 
T300/5208 graphite-epoxy 
[45/0/-45/90] g laminate 
Rectangular 90 mm x 150 mm 


Impact 

Rebound 

Coefficient of restitution, e 

Drop 

height, m 

Maximum 
velocity, m/s 

Rebound 
height, m 

Initial 
velocity, m/s 

Test 

Computed 

0.30 

2.43 

0.21 

-2.03 

0.836 

0.8355 

.60 

3.43 

.35 

-2.62 

.764 

.7825 

.90 

4.20 

.57 

-3.34 

.796 

.7464 

1.20 

4.85 

.67 

-3.62 

.747 

.7184 

1.50 

5.43 

.87 

-4.14 

.762 

.6951 

1.83 

5.99 

.92 

-4.25 

.709 

.6752 


TABLE II 

RESULTS OF CANTILEVER IMPACT TESTS 

Impactor : 1-inch diameter steel ball 

Target : T300/5208 graphite-epoxy 

[45/0/-45/90] g laminate 

Impact area: Circular 90 mm and 50 mm diameter 



D = 

90 ram 



D = 

50 mm 


V 

v r 


e 

V 

V_ 


e 

° 




0 

r 



m/s 

m/s 

Test 

Computed 

m/s 

m/s 

Test 

Computed 

0.48 

-0.460 

0.960 

0.962 

0.40 

-0.400 

1.0 

0.968 

.74 

-.705 

.952 

.943 

.725 

-.674 

.930 

.944 

1.06 

-.943 

.890 

.920 

.930 

-.825 

.887 

.929 

1.22 

-1.146 

.939 

.909 

1.200 


.950 

.911 

1.52 

-1.368 

.900 

.890 

1.510 

■ 

.876 

.891 

1.98 

-1.665 

.841 

.861 

2.600 

-2.189 

.842 

.826 

' 2.60 

-1.926 

.741 

.325 

3.150 

-2.598 

.825 

.796 I 

3.39 

-2.779 

.820 

.784 

3.690 

-2.911 

.789 

.770 

3.60 

-2.898 

.805 

.774 

4.210 

-3.200 

.760 

.750 

3.92 

-3.136 

.800 

.761 

5.300 


.733 

.700 

4.98 

-3.790 

.761 

.713 





5.47 

-3.993 

.730 

.693 





6.12 

-4.437 

.725 

.673 






V Q , Impact velocity 

V r , Rebound velocity 

e. Coefficient of restitution 
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Caption 


I. . Idealized dynamic behavior of the laminate. 

2. Variation of impact affected zone size with impact velocity. 

3. Effect of shear modulus on the size of impact affected zone. 

4. Effect of interlaminar shear strength on impact affected zone size. 

5. Effect of composite's mass density on impact affected zone size. 

6. Variation of coefficient of restitution (C.O.R.) with impact velocity. 

7. Effect of shear modulus on coefficient of restitution. 

8. Effect of interlaminar shear strength on coefficient of restitution. 

9. Effect of composites mass density on coefficient of restitution. 

10. Test set-up for free ball drop impact tests. 

II. Test set-up for cantilever type instrumented impact tests. 

12. Variation of impact force, velocity and deflection with time. 

13. Comparison of predicted and test values of coefficient of restitution. 
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Effect of composite's mass density on impact affected zone size. 
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7. Effect of shear modulus on coefficient of restitution. 
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10. Test set-up for free ball drop impact tests. 
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Variation of impact force, velocity and deflection with time. 
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APPENDIX B 


Relationship Between Phase Difference and Coefficient of 
Restitution During Low Velocity Foreign Object Transverse 
Damage of Composite Plates 
K. M. Lai* 

NASA Langley Research Center 
Hampton, Virginia 


Abstract: 

This work discusses a model to correlate the coefficient of restitution, 
of low velocity transverse impacts of graphite-epoxy laminates, with the 
residual deformation or central deflection at the end of the impact event. 

It is assumed that the energy absorbed by the target can be calibrated in 
terms of residual deflection, and thereby in terms of phase difference 
between the occurrence of impact force and central deflection to their 
maximas. Analysis was modeled on the basis of the experience from impact 
tests. 

Predictions are compared with the test results of impacted circular 
and flat, plates. Experimentally measured, values of coefficient of restitu- 
tion and phase difference agreed well with the predicted relationship 
between them. 


*Also, Research Associate Professor of Mechanical Engineering and Mechanics at Old 
Dominion University, Norfolk, Virginia. 
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During the last decade, significant progress has been made in develop- 
ing advanced composites possessing high strength, high modulus, and low 
density, and in understanding their behavior under certain types of loading. 
As a result, composites such as graphite-epoxy and boron epoxy have been 
successfully employed as structural materials in aircraft, missiles, and 
space vehicles, and have satisfactorily demonstrated their performance 
through extensive ground testing and in flight. 

Despite the tremendous advantages that advanced composites have over 
metals in applications requiring high strength, high stiffness, and low 
weight, in applications where impact by foreign objects is a design con- 
sideration, the advantages inherent in the composites are overshadowed by 
their poor response to impact loading. Aircrafts, while in service, may 
be subjected to impacts from hails and bird strikes in the air, runway 
debris, and even ground service equipment. Such impacts can reduce the 
strength of the structure. 

Literature survey [l] of low velocity foreign object damage to com- 
posites, in particular the thin plates of graphite-epoxy, has revealed that 
the following mechanisms may be expected to operate during the impact event: 

(1) Pullout 

(2) Debonding 

(3) Post debonding friction work 

(4) Surface energy 

It should be emphasized that it is not necessary that all of the above 
four mechanisms to operate simultaneously for a given fiber-matrix system. 

As a result of the above mechanisms, the energy absorbed by the target 
during impact can be related by the coefficient of restitution [2] such that 


abs 


1(1 - e 2 ) 


( 1 ) 
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where I is impact energy and e is coefficient of restitution of the 
target. In this work, an attempt is made to correlate the factor e with 
the phase difference between the maximum impact force and central deflection 
as seen in the tests. 

Tests 

In order to understand the impact mechanism of composite plates, few 
impact tests simulating the low velocity impact conditions were performed. 
Specimens were of circular plates of 50 mm and 100 mm diameter with clamped 
boundary. 8-ply and 16-ply plates were tested. Ply orientation, 
respectively, for 8 -ply and 16-ply plates were [45,0,-45,90] g and 
[452»02»-452,902] g . Impact was given by an instrumented impactor fixed on 
the free end of the cantilever. Fig. 1. During an impact event, the elec- 
trical signals from the accelerometer (clamped on the top of the impactor) 
were first processed through a digital data acquisition system, and then 
through a computer and placed onto cassette tapes to allow direct 
acceleration-time data plotting. Double integration of the test data 

\ 

provided respectively the velocity-time and displacement data. Typical 
variation of impact force and displacement with time is shown in Fig. 2. 

It can be seen that the variation of force looks like sinusoidal, it is zero 
at the beginning and at the end of the impact event, and that it is maximum 
at approximately half of the impact duration. The variation of central 
displacement is also sinusoidal with a phase lag with the impact force. 

The deflection was zero at the beginning of the impact, reaches to its 
maximum after the force reaches its maximum value, and there is some 
residual deflection at the end of the impact event. Amount of residual 
deflection was increasing with severeness of the impact. This behavior is 
described that the applied load and resultant deflection are out of phase. 
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As a result of this phase difference, there xs internal friction absorbing 
energy during the impact event [3]. Obviously, there may be external 
frictional effects such as air damping, but in present discussions we shall 
assume these to be negligible and that all damping is internal. 

Based on the experience from actual impact test behavior, the variation 
of impact force and central deflection with time during an impact event are 
modeled in the next section. 

Analysis 

The amount of damping or the magnitude of internal friction may be 
measured in a number of ways, one such being the phase angle <f> . 

Consider the load applied during the impact event has the form 

P(t) = P q sin £ t (2) 

where P(t) is the load at any time t, T is impact duration, and P q is 
maximum load during impact. The load will produce a periodic displacement, 
D(t) at the point of impact, such that 

D(t) = D q sin ^ ~ ^ ]t (3) 

where D(t) is the displacement at time t, D q is maximum deflection, 
and <J> is a parameter for phase difference angle. Putting t = T in 
equation (3) gives the residual deflection, D r , at the end of the impact 
event : 


D = D sin <j)T 
r o T 

4>T = sin _1 (D r /D o ) 


or 


(4) 
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Here <J)T is the measure of phase difference, between the maximum load and 
maximum deflection, in radians. 

The work done or the energy absorbed by the target during impact event 
is 


I 


abs 



P(t) d (D (t) ) dt 


(5) 


Substitution of P(t) and D(t), respectively, from equations (2) and (3) 
and integrating from zero to T gives 


I 


abs 


P D (1 - z) 

O 0 V „ _ X 

z ( 2 z ) ~ (1 ■ COS Uz) 


( 6 ) 


where z = ^T/lT. The time, t Q , when the central deflection is maximum can 
be obtained from equation (3) 


to-ir/a-z) 


(7) 


The maximum energy, I, of the target is when the deflection is 

maximum. This is obtained by integrating the equation (5) from time 

t = zero to t , and is 
o’ 


The ratio of I , 
abs 

of restitution, e, 
following equation: 


(1 - z) 


, / TTZ 

\ 

"rr(2 - z)\ 

1 C0S \2(1 - z) 

J 1 cosl 

.2(1 - z )/ 

z 

(2 

- z) J 


and I gives the relative damping. The coefficient 

can be written in terms of the ratio I , /I by the 

abs 


( 8 ) 



e = 


(9) 
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Substitution of I an< * I, respectively, from equations (6) and (8) into 
equation (9) gives 



( 10 ) 


The values of coefficient of restitution, e, were computed from 
equation (10) and are plotted in Fig. 3 as a function of phase difference. 
It is interesting to mention that the variation of coefficient of restitu- 
tion with phase difference is independent of time duration of impact. 


Results and Discussions 

The impact data obtained from tests are given in Table I. Data are 
plotted on Fig. 3 to compare the validity of the prediction. The agreement 
between theory and experiment was found better than could reasonably be 
expected, in view of very simple nature of the theory. 

It is important to mention that the predicted relationship between 
coefficient of restitution and phase difference agrees well with the results 
of flat surface impact, as well as of clamped plate impact tests. This 
suggests that the energy absorbed by plates could be predicted by impact 
testing the flat surface specimens. 
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CAPTIONS 


Fig. 1. Lay-out of the cantilever- type instrumented impactor. 

Fig. 2. A typical response of acceleration and deflection obtained from 
accelerometer . 

Fig. 3. Comparison of predicted relationship between phase difference and 
coefficient of restitution with the test data. 
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Fig. 3. Comparison of predicted relationship between phase difference and 
coefficient of restitution with the test data. 
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